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Abstract 

Self-doping in graphene has been studied by examining single-layer epitaxially-grown graphene 
samples with differing characteristic lateral terrace widths. Low energy electron microscopy was 
used to gain real-space information about the graphene surface morphology, which was compared 
with data obtained by angle-resolved photoemission spectroscopy to study the effect of the mono- 
layer graphene terrace width on the low energy dispersions. By altering the graphene terrace width 
we report significant changes in the electronic structure and quasiparticle relaxation time of the 
material, in addition to a terrace width-dependent doping effect. 

PACS numbers: 
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Single-layer graphene has been the subject of much research in recent years jl| due to 
the remarkable physics it contains and its potential for device applications. The electronic 
properties of graphene sheets are extremely sensitive to their boundary conditions: chang- 
ing these boundary conditions can create and alter semiconducting band-gaps [^j], can lead 
to half-met allicity, 3| and precise control over graphene edge geometries may allow for the 
creation of complex device components J4]. Such control is even more important in light 
of the fact that edge geometries and structural disorder can lead to the formation of im- 
purity bands near the Fermi level, [5] effectively hole-doping the system and substantially 
altering the transport properties of graphene. Although the effect of this self-doping in 
ribbons has been analysed theoretically, fl so far no systematic cental study e xis t, 
Transport studies of graphene nanoribbons have been performed, [2|, [6|] but the p-doping 
may be attributed to the patterning process and atmospheric contaminants. Similarly, al- 
though evidence of self doping can be found in Raman measurements, 

m 

the doping 

levels vary greatly among these samples even in the absence of measured disorder. The pri- 
mary difficulty is in the comparison of the microscopic features of disorder to spectroscopy 
data: the size of the probing beam is often too large to study the position-dependence of 
the self-doping as a function of the distance from a sample edge. Instead, controlling the 
characteristic length scale of the sample itself can be a more practical solution. This can be 
achieved by growing epitaxial graphene on 6H-SiC, as it forms a tangle of broad nanoribbon- 
like monolayer graphene terraces 10|, 111] , whose size can be directly monitored by low energy 
electron microscopy (LEEM) 12| during the growth process. The effect of disorder and/or 
self-doping on the electronic structure of graphene can then be directly studied on the same 
samples by angle-resolved photoemission spectroscopy (ARPES) 13j. 

In this letter we present a systematic study of the self-doping on epitaxially-grown 
graphene through the edges of the graphene terraces. Several epitaxially-grown graphene 
samples with average monolayer terrace widths ranging from 50nm to 180nm have been 
studied by LEEM and ARPES. As the terrace widths increase we observe an increase in 
charge carrier density, resulting in a rigid shift of the low binding-energy bandstructure and 
a constant Fermi velocity, and a corresponding decrease of the quasiparticle lifetime due to 
the increase in disorder. These results provide direct experimental evidence of self-doping 
effects in graphene and illuminate the role of finite graphene ribbon widths on its electronic 
structure. 
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Samples were grown on the Si-terminated face of 6H-SiC substrate as previously reported 
[ijj]. LEEM was performed at the National Center for Electron Microscopy at Lawrence 
Berkeley National Laboratory lo[ 11]. High resolution ARPES data were taken at BL12.0.1 
of the Advanced Light Source at a temperature of 15K with 50eV photons after annealling 
samples to 1000C. 

Figures l(a-c) show LEEM images of three characteristic ep it axially- grown graphene 



samples. 
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11| The width of the single-layer graphene terraces has been extracted using a 
standard method 3], by first identifying outlines of the single-layer terraces, then drawing 
straight lines in the images, and finally measuring the linear widths of the terraces crossed 
by the lines. The terrace size increases from 50nm to 180nm from panel a to panel c. 

In figure 2, the raw ARPES intensity maps (lower panels) and area-normalized momentum 
distribution curves (MDCs), intensity profiles at constant energy (upper panels), are shown 
for five samples, S1-S5, ordered by increasing terrace width. Data are taken along a high 
symmetry direction through the K point (see horizontal black line in the inset of panel a). 
As the terrace widths increase the following qualitative trends are observed (quantitative 
discussion follows in figures 3 and 4): 1) The separation between the two branches that 
disperse towards E^ (dotted line) increases (compare the separation between the MDCs 
peaks in the upper panel); 2) The intensity maps become sharper at the Fermi level and 
higher binding energy (compare width of the MDCs peak in the upper panel); 3) The region 
of vertical intensity separating the conduction band from the valence band, namely the 



gap : 



16] decreases in agreement with a previous report. [17| It can also be seen in panels d 



and e that additional bands appear faintly (see bottom of the conduction band and a replica 
of the valence band) due to an increase in the ratio of the bilayer to the monolayer content 
as its coverage increases. 

A quantitative analysis of the Fermi level MDCs is presented in Figure 3. For ease of 
comparison, the MDCs curves shown in the previous figure have been offset and normalized 
by the area underneath, after removing a constant background (see panel a). The direct 
comparison clearly shows a shift of the MDC peak positions toward higher momentum 
and a decrease of the MDC widths with increasing terrace size, pointing to a change of the 
doping level and scattering rate, respectively. The charge density (panel b) and quasiparticle 
scattering rate (panel c) as a function of the monolayer graphene terrace widths are shown. 
The charge density, a quantity directly related to the separation of the MDCs peak at 
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Ef, is calculated by assuming the electron pockets of graphene to be circles in k-spacefis| 
with diameters given by the separations of the dispersions at the Fermi level, and invoking 
Luttinger's theorem to extract the charge carrier density from the Fermi surface area. Clearly 
the size of the electron pocket decreases, or the hole doping increases, as the monolayer 
graphene terrace width increases from 50nm to 180nm, in line with the shift of the MDC 
peak position of panel a. This is in agreement with the prediction that self-doping varies 
inversely with the graphene ribbon width[s|. More specifically, the data can be well described 
by the function a = CiW -1 + C2 (gray dashed line in panel b), where a is the charge density, 
W is the terrace width, and c\ and C2 are constants. We find that the constant Ci=-2xl0 7 
cm -1 is very close to the approximate theoretical value of -5xl0 6 cm -1 (using t'=0.1eV, and 
a dielectric constant of 10 in Eq. 161 of ref jf]]). In the limit of an infinite epitaxially-grown 
graphene plane, one would expect an electron-doping of approximately C2 = 1.17xl0 13 cm -2 
due to the interaction with the substrate. 

Despite the overall agreement between theory and experiment, the theoretical prediction 
for the change in doping as the ribbon width changes from 34nm to 170nm is 6.5 xlO 11 
electrons/cm 2 , which is a factor of ~4 smaller than the experimentally observed doping for 
similar widths. There are several reasons why this might be expected. First, the method 
used for extracting the terrace widths from LEEM images is likely an overestimate of the 
actual ribbon width, since an arbitrary line drawn across a LEEM image will not cross the 
ribbon edges perpendicularly, and because the resolution of the LEEM is finite. Second, the 
experimental geometry is different from the theoretical one: jagged edges might permit more 
self-doping than straight ones, and higher densities of other types of lattice defects may be 
associated with real ribbon edges. Finally, the details of the self-doping of epitaxially-grown 
graphene might be different from the free-standing case. 

Figure 3(c) shows the Fermi level MDC width, a quantity directly related to the inverse 
quasiparticle relaxation time, T ~ 1/r, as a function of the monolayer graphene terraces 
width. As the terrace widths increase we observe a clear decrease of the MDC width, 
or increase of the quasiparticle scattering rate, in line with panel (a) (compare horizontal 
arrows). If the quasiparticle scattering rate were dominated by the terrace edges, the mean 
free path, and therefore the transport lifetime, would vary linearly with the terrace width, 
although it should be noted that the single-particle relaxation time probed by ARPES may 
not equal the transport lifetime [191] . Therefore, as a guide to the eye, a dotted line is plotted 
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of the form FWHM = c^W' 1 + c 4 , where FWHM is the MDC width, W is the terrace 
width, and C3 and C4 are constants equal to 7.5 and 0.0125A -1 , respectively. This change in 
lifetime reflects the decrease of the quasiparticle mean free path with decreasing monolayer 
graphene terrace width and provides further confirmation of the observation of self-doping 
in graphene. 

The effect of finite terrace size on the conical dispersion is shown in figure 4. The energy 
vs momentum dispersion relation are obtained by fitting the MDCs spectra over the full 
energy range, using the standard Lorentzian fitting procedure. [3] Dispersions for binding 
energies lower than the Dirac point (\E — Ep\ < \E — Ed\) are displayed in panel a, while 
dispersions at higher binding energy (\E — Ep\ > \E — E£>\) are displayed in panel b. With 
increasing terrace width the data show a clear rigid shift of the dispersion away from the 
K point at low binding energy (in line with the observed shift of the MDC peaks at kp 
in panel a) and a rigid shift of the dispersions toward the K point at high binding energy. 
These shifts, together with the nearly constant Fermi velocity (panel c) extracted from the 
low energy dispersions Vp=|^j|, point to a rigid shift of the conduction band in energy, and 
consequently an increase of the chemical potential, as the terrace widths increase. For a 
schematic summary see inset of panel a. Interestingly, although on a qualitative level the 
shift observed in the high energy dispersion of panel b is also in agreement with this picture, 
on a quantitative level the magnitude of the shift is smaller by roughly a factor of 2. This 
apparent discrepancy can bepartially accounted for by the reported change in the size of 



the gap at the Dirac point 



171 ]. but further investigation is required to exactly determine 



the mechanism behind this change. 

In conclusion, we have manipulated the electronic structure of monolayer graphene by 
adjusting the widths of the monolayer terraces during the growth process of epitaxially 
grown graphene. We found that altering the graphene terrace width has a significant effect 
on the electronic structure of the material, even for samples with length scales on the order of 
lOOnm. These finite-size effects will play a role in electronics applications as graphene-based 
components of various sizes and geometries are fabricated in the future. 
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Figure captions: 



Fig. 1: (Color online) (a)-(c) 2/xm x 2/im LEEM images of three samples, with red dots to 
denote the typical terrace widths used in this paper. The three shades in these images, white, 
grey, and black, are buffer layer, monolayer graphene, and bilayer graphene, respectively. 
Inset of (a): Reflectivity as a function of incident kinetic energy for the buffer layer and first 
monolayer through fourth monolayer (bottom to top). The number of dips in the reflectivity 
curves is equal to the number of graphene layers. The dotted line denotes 4.2eV electron 
kinetic energy, which is used in (a)-(c). 

Fig. 2: (Color online) (a-e) ARPES spectra of samples with increasing monolayer graphene 
terrace width, SI through S5. Inset of (a) shows the Brillouin zone of graphene. The 
horizontal line through the K point shows the measurement geometry in k-space. The large 
and small pink left-right arrows give the peak separation and width of the Fermi level MDC 
of SI. 

Fig. 3: (Color online) (a) The Fermi level MDCs with incremental offset, (b) The electronic 
charge carrier density as a function of terrace width, (c) Fermi level MDC widths, extracted 
from panel a. 

Fig. 4: (Color online) Low (a) and high (b) binding energy dispersions as a function of 
terraces widths. Note that the small curvature at E^ in panel a is due to the finite energy 
resolution of the experiment, (c) Fermi velocities vs terrace widths. Inset of (a, b) shows a 
cartoon of the graphene band structure as a function of terraces width. The boxed regions 
in the insets denote the regions of the dispersions of Figs. 4(a) and 4(b). 
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FIG. 1: 
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